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Experimental results showing huge negative differential conductance in gold-hydrogen molecular 
nanojunctions are presented. The results are analyzed in terms of two-level system (TLS) models: it 
is shown that a simple TLS model cannot produce peaklike structures in the differential conductance 
curves, whereas an asymmetrically coupled TLS model gives perfect fit to the data. Our analysis 
implies that the excitation of a bound molecule to a large number of energetically similar loosely 
bound states is responsible for the peaklike structures. Recent experimental studies showing related 
features are discussed within the framework of our model. 
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I. INTRODUCTION 

The study of molecular nanojunctions built from 
simple molecules has attracted wide interest in recent 
years ^ Contrary to more complex molecular electron- 
ics structures, the behaviour of some simple molecules 
bridging atomic-sized metallic junctions can already be 
understood in great detail, including the number of 
conductance channel analysis with conductance fluctu- 
ation and shot-noise measurements^^ the identifica- 
tion of various vibrational modes with point-contact 
spectroscopy^ and the quantitative predictive power of 
computer simulations^ The above methods were success- 
fully used to describe platinum-hydrogen junctions: it 
was shown that a molecular hydrogen bridge with a sin- 
gle, perfectly transmitting channel is formed between the 
platinum electrodes^ 

Point-contact spectroscopy turned out to be an espe- 
cially useful tool in the study of molecular nanojunctions, 
as a fingerprint of the molecular vibrational modes can be 
given by simply identifying the small stcplike signals in 
the dI/dV(V) curves of the junction^ However, the de- 
tection of the small vibrational signals is difficult for junc- 
tions with partially transmitting conductance channels, 
where quantum interference (QI) fluctuations may give 
an order of a magnitude larger signal. 6 Surprisingly, in- 
stead of observing small steplike vibrational signals upon 
the background of QI fluctuations, molecular nanojunc- 
tions frequently show peaklike structures in the differen- 
tial conductance curves with amplitudes comparable to 
or even much larger than the QI fluctuations. The peaks 
can be either positive or negative and their amplitude can 
be as small as a few percents, but - as we demonstrate in 
this manuscript - the peak-height can be several orders 
of magnitude larger showing a huge negative differential 
conductance (NDC) phenomenon (Fig.[T]). 

As a general feature of the phenomenon it can be stated 
that at low bias the conductance starts from a constant 
value, at a certain threshold voltage the peaklike struc- 
ture is observed, and at higher bias the conductance sat- 
urates at a constant value again. The high-bias conduc- 
tance plateau can be either higher or lower than the low- 
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FIG. 1: Huge negative differential conductance (upper panel) 
and the corresponding I — V curve (lower panel) in gold- 
hydrogen nanojunctions. The solid black lines show the theo- 
retical fits with the asymmetrically coupled TLS model. The 
parametrs are: E = 54meV, T — 5.1 K, <xo = 0.45 Go, 
(Too = 0.12 Go, N = 60 andW = 0. 



bias plateau. In the first case the peak at the transition 
energy is positive, whereas in the second case it is neg- 
ative. In other words, the junction shows linear I — V 
characteristics both at low and high biases but with a 
different slope, and at the threshold energy a sharp tran- 
sition occurs between the two slopes. The switching of 
the conductance between two discrete levels implies the 
description of the phenomenon with a two-level system 
(TLS) model, but as we demonstrate in this paper a sim- 
ple TLS model only gives steps in the dl/dV at the exci- 
tation energy of the TLS, and it cannot account for sharp 
peaks. 

The above phenomenon is not a unique feature of a 
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special atomic configuration, but it frequently occurs in 
a wide conductance range and it appears almost in all the 
studied molecules and contact materials. In our exper- 
iments on gold-hydrogen molecular junctions^ we have 
recognized huge negative differential conductance curves 
in the conductance regime of 0.1 — 0.7 Go showing a char- 
acteristic threshold energy of 30— 100 meV. We have also 
observed similar phenomenon in niobium-hydrogen nano- 
junctions. In parallel the group of Jan van Ruitenbeek 
- focusing on the conductance range close to the quan- 
tum conductance unit - has demonstrated the occurrence 
of smaller peaklike structures^ In their study the peaks 
appeared with H 2 , D 2 , O2, C2H2, CO, H 2 and ben- 
zene as molecules, and Pt, Au, Ag and Ni as contact 
electrodes. An scanning tunneling microscopy (STM) 
study of a hydrogen-covered Cu surface revealed both 
sharp gaplike positive peaks and large negative differen- 
tial conductance peaks in the tunneling regime^ The lat- 
ter study has also demonstrated a clear telegraph fluctu- 
ation between the two levels on the millisecond timescale. 
An earlier STM study demonstrated smaller negative dif- 
ferential conductance peaks due to the conformational 
change of pyrrolidine molecule on copper surfaced 

The above results imply the presence of a general phys- 
ical phenomenon of molecular nanojunctions resulting in 
a similar feature under a wide range of experimental con- 
ditions. In Ref. 8 it was shown that in certain systems the 
peaklike structures are related to the vibrational modes 
of the molecular junctions and the observed features were 
successfully described in terms of a vibrationally medi- 
ated two-level transition model. This model, however, 
cannot describe the large negative differential conduc- 
tance curves in our studies, and the large peaklike struc- 
tures in the STM experiments^ 

In the following, we present a detailed analysis of the 
observations in terms of two-level system models. We 
demonstrate the failure of a simple TLS model and the 
necessary ingredients for producing peaklike structures 
instead of simple conductance steps. We propose a model 
based on an asymmetrically coupled two-level system, 
which can describe all the above observations. The com- 
parison of the model with experimental data implies a 
huge asymmetry in the coupling strength, which can 
be explained by exciting a strongly bound molecule to 
a large number of energetically similar loosely bound 
states. In our opinion, the proposed model is appropri- 
ate for describing the appearance of peaklike structures 
(or even negative differential conductance) in the dl/dV 
curves of molecular nanojunctions under a wide range of 
experimental conditions. 



II. THE FAILURE OF A SIMPLE TWO-LEVEL 
SYSTEM MODEL 



ences therein). In the following, based on the results for 
a general point-contact geometry^ we shortly give an 
overview of the scattering process on a two-level system 
located near the center of an atomic-sized nanojunction. 
The TLS is considered as a double well potential, where 
the two states in the two potential wells have an energy 
difference A (see Fig. ^ . The two states are coupled by 
tunneling across the barrier with a coupling energy T. 
The coupling between the wells causes a hybridization of 
the two states, resulting in two energy eigenvalues with 
a splitting of E — %/A 2 + T 2 . It is considered that in the 
lower state of the TLS the contact has a conductance ao, 
whereas in the upper state the conductance is o\. The 
occupation number of the two states are denoted by no 
and n\ = 1 — no- These are time-averaged occupation 
numbers, the system is jumping between the two states 
showing a telegraph fluctuation. On a time-scale much 
longer than that of the TLS the fluctuation is averaged 
out, and the I — V characteristic is determined by the 
voltage dependent occupation numbers: 



I(V) = (<7 7io + °in>i)V, 
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The occupation of the upper state can be calculated from 
the rate equation: 
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where the transition probabilities are determined by 
Fermi's golden rule: 

JW = p|no7 / def(e, eV)(l - f(e - E, eV)), (4) 
Pi^o = p F n xl f de/(e, eV)(l - f(e + E, eV)). (5) 



Here an electron from an occupied state with energy e 
scatters on the TLS to an unoccupied state with energy 
e ± E. The nonequilibrium distribution function of the 
electrons is denoted by f(e,eV), whereas pp stands for 
the density of states at the Fermi energy. The transi- 
tion matrix element from an initial electron state and the 
ground state of the TLS to a final electron state and the 
excited state of the TLS is considered as a constant cou- 
pling strength: 7 = 2%/h- 0\H e - T Ls\f, 1) | 2 - Assum- 
ing that the TLS is situated at the middle of the contact, 
where the half of the electrons is coming from the left 
and the half from the right electrode, the nonequilibrium 
distribution function can be approximated as: 



f(e,eV) = 



(6) 



The scattering on two-level systems in mesoscopic 
point-contacts has been widely studied after the discov- 
ery of point-contact spectroscopy (see Ref. [Ill a- n d refer- 



where /o'(e) and /^(e) = /q (e— eV) are the equilibrium 
Fermi functions of the left and right electrodes, for which 
the chemical potential is shifted by the applied voltage. 
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By inserting the distribution function to Eq. [4][5] and us- 
ing the formula: 



/ ( e )(l-/ (e-a))d6 = |(coth^ 



the rate equation can be written as: 

dm 



i) 



dt 



= n Q v -nxVi, 



(7) 
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where vq and v\ are the inverse relaxation times of the 
lower and upper state of the TLS, 
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The steady state solution for the occupation number of 
the upper state is: 



Hi = 1 — n = 



(11) 
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FIG. 2: Results in the zero temperature limit of a symmetri- 
cally coupled TLS model. The upper panel shows the voltage 
dependence of the occupation numbers. The characteristic en- 
ergy describing the growth of ni above the excitation energy is 
denoted by SE. The middle panel shows the I — V curve, the 
low and high bias slopes (<tq - V, a x - V) are illustrated by dot- 
ted lines. The lower panel shows the differential conductance 
exhibiting a jump of Ao\e at the excitation energy and an 
overall change of Acr,^ between the zero and high bias limits. 
On the right side the two-level system model is demonstrated. 



In the zero temperature limit the inverse relaxation times 
and the occupation number take the following simple 
forms: 
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Fig. [2] shows the evolution of the occupation numbers, 
the I — V curve and the differential conductance curve 
in the zero temperature limit. At eV ^S> E both states 
are equally occupied with n = 1/2, but the transi- 
tion towards this is very slow. A characteristic energy 
scale, SE for the variation of the occupation numbers 
can be defined by extrapolating the linear growth of n\ 
at eV = E + to the saturation value (see the upper panel 
in the figure). Due to the small slope of ni(V) this char- 
acteristic energy scale is larger than the excitation energy, 
precisely: 8E = 2E. This slow change is reflected by a 
smooth variation of the I — V curve. The differential 
conductance curve shows a steplike change at eV = E, 
and then saturates to CTqc = (ao + <Ji)/2. The size of 
the step is smaller than the overall change of the conduc- 
tance (Ao\e = AfJoo/2), thus only a step and no peak is 
observed in the dl/dV curve. 



In the above considerations the calculation of the cur- 
rent (Eq. [T|) includes only elastic scattering of the elec- 
trons on the two states of the TLS with different corre- 
sponding conductances, and the inelastic scattering pro- 
cesses are just setting the voltage dependence of the oc- 
cupation numbers. However, the inelastic scattering of 
the electrons in principle gives direct contribution to the 
current: above the excitation energy the backscattering 
on the contact is enhanced due to the possibility for in- 
elastic scattering. More precisely the inelastic current 
correction can be written as: 

SI m = e ■ [(-P l+1 o-l_,i + Pl + ,o^r + ,i (15) 
— -Pr_,o->l_,i + -Pr_,o^r + ,i) + (o«i)]' 

Here L + and R_ are incoming states to the contact on the 
left and right side, while R + and L_ are outgoing states, 
and e.g. -Pl + ,o-*l_,i is the probability that an electron 
from an incoming state on the left side excites the TLS 
and scatters to an outgoing state on the left side. The 
incoming states have the distribution function of the cor- 
responding reservoir, while the distribution functions of 
the outgoing states are a mixture of the left and right 
Fermi functions with the transmission probability of the 
contact, t: 

The sign of the different current terms are determined by 
the direction of the outgoing states. After evaluating the 
energy integrals with the appropriate distribution func- 
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tions one obtains: 
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where the first and second term correspond to the exci- 
tation and relaxation of the TLS, and: 
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The inelastic correction also causes a steplike change 
of the conductance, which has a magnitude of SG m « 
e 2 p|,7(l — 2r)/4. For a contact with large transmission 
the conductance decreases at the excitation energy, while 
for a tunnel junction it increases, with a transition be- 
tween the two cases at r = 1/2. We note that our simple 
model gives the same result for the transition between 
point-contact spectroscopy and inelastic electron tunnel- 
ing spectroscopy at r = 1/2 as more detailed calculations 
for vibrational spectroscopy signals i 12 i 13 

The contribution of the inelastic process can easily be 
estimated. According to Eq. [13] the relaxation time of 
the TLS is r TLS ss {ppjE) -1 , thus the inelastic correc- 
tion to the conductance is \SG m \ < e 2 /(r T LS • E). It 
means that a relaxation time of 1 /isec corresponds to an 
inelastic correction smaller than 10 -7 Go- In other words 
a correction of 0.1 Go would correspond to a TLS with 
a sub-picosecond relaxation time. These are unphysical 
numbers, especially if the telegraph fluctuation can be 
resolved, thus in the following the inelastic correction is 
neglected. 

As a conclusion a simple TLS model cannot produce 
pcaklike structures in the differential conductance curve. 
Regardless of the fine details of the model, for any TLS 
a characteristic energy 6E > E is required for the tran- 
sition, whereas for the observation of sharp peaks in the 
differential conductance 5E <C E is desired. As an other 
consequence of a simple TLS model, the high voltage 
conductance is limited by CToo > cto/2 due to the equal 
population of the states at high bias. This contradicts 
the observation of Fig. [T] where <to/5. For such 

large change of the conductance a population inversion, 
ni(oo) ^> no(oo) is required. 



III. TWO-LEVEL SYSTEM MODEL WITH 
ASYMMETRIC COUPLING 

Both a population inversion and a sharp transition of 
the occupation numbers can be introduced by inserting 
an asymmetric coupling constant to the model, that is the 
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FIG. 3: Results in the zero temperature limit of an asymmet- 
rically coupled TLS model. The upper panel shows the volt- 
age dependence of the occupation numbers, the middle panel 
shows the I — V curve, and the lower panel shows the differen- 
tial conductance. On the right side the two-level system with 
a degenerated upper level is demonstrated. 



coupling of the lower state of the TLS to the electrons in 
Eq. [4] is 70, whereas the coupling to the upper state in 
Eq. [5] is 71 <C Jo. The asymmetry parameter is defined 
as N — 7o/yi- With this modification, once the TLS 
can be excited, it cannot relax back easily, thus a sharp 
transition occurs. 

The original definition of the coupling constant 
does not allow any asymmetry, as 70 = 2n/h ■ 

|(/,l|# e _TLsM)| 2 = 27T/R- |(z,0|ff e _TL5|/,l)| 2 = 71 

due to the hermicity of the Hamilton operator. An asym- 
metry arises, however, if the phase spaces of the two 
levels are different. For instance if the ground state is 
well-defined, but the upper state is not a single level, but 
N energetically equivalent states (Fig. [3]) , then the effec- 
tive coupling constant contains a summation for the final 
states resulting in 70 = A/71 . 

With this modification the occupation numbers in the 
T = limit are defined by: 
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The corresponding occupation numbers, I—V and dl/dV 
curves are plotted in Fig. [3l The population inversion 
is obvious with ni(oo)/nn(oo) = N. The characteris- 
tic energy of the variation of the occupation numbers is 
6E = AE/ (N + 1). The I—V curve already shows a 
sharp transition from the initial to the final slope, and 
the differential conductance shows a sharp peak if the 
asymmetry parameter is high enough. The appearance 
of the dl/dV peak is determined more precisely by cal- 
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With an asymmetry parameter N > 3 a peak is observed 
in the differential conductance, whereas for smaller asym- 
metry only a steplike change arises. For a given differ- 
ence in the conductances, (a± — <r ) both the width and 
the height of the conductance peak are determined by 
the asymmetry parameter, i.e. the peak width and the 
peak hight are not independent parameters. Depend- 
ing on the sign of (a% — <jq), the peak can be either 
positive or negative. The conductance at large bias is 
Coo = {N<T\ + gq)/{N + 1), thus it can be arbitrarily 
small compared to the zero-bias conductance. 

At finite temperature the dl/dV curve is calculated 
by inserting Eqs. [TOl HU into Eq. [5] using asymmetric 
coupling constants, 70, 7i- The finite temperature causes 
a smearing of the curves, but the transition from steplike 
to peaklike structure is similarly observed at N = 3. 

A more realistic generalization of the model can be 
given by assuming a distribution of the energy levels at 
the excited state with a density of states pi(E), for which 
the standard deviation (W) is kept much smaller than the 
mean value (Eq). The asymmetry parameter is defined 
by the normalization: N = j pi(E)dE. With this mod- 
ification the occupation of the excited states is energy 
dependent: rii(E), but the probability that any of the 
excited states is occupied is given by a single occupation 
number, n\ = J n\{E)dE. The calculation of the tran- 
sition probability Po— >i includes the integration of the 
excitation rate in Eq. [9] with the energy distribution: 



u (eV,E ,W) = / Pl (E)u (eV,E)dE. 



(22) 



The precise calculation of the relaxation probability 
(Pi_>o) would require the knowledge of the energy depen- 
dent occupation number, m(E). However, in the narrow 
neighborhood of the excitation energy, where the peak is 
observed the voltage dependence of the relaxation rate 
can be neglected beside the constant value of w AE (see 
Eq.[T3]), thus the original formula for v x (Eq.QI)]) is a good 
approximation for the finite distribution of the levels as 
well. Our analysis shows, that the relaxation rate can 
even be replaced by a voltage independent spontaneous 
relaxation rate without causing significant changes in the 
results. 

As demonstrated in the following a TLS model with 
an asymmetry parameter, N, and a narrow width of the 
upper states, W, shows perfect agreement with the ex- 
perimental observations. 
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FIG. 4: Panel (a) shows a conductance trace on which NDC 
curves were observed when the electrodes were approaching 
each other after rupture. In this region the differential con- 
ductance measured at low DC bias (ao, Vdc = OmV) and 
high bias ((Too, Vdc = 120 mV) exhibit a large splitting. The 
differential conductance was measured with an AC modulation 
of 500 fiV. Panel (b) shows conductance histograms of Au-Hi 
junctions measured at a DC bias of 20 mV and 120 mV, re- 
spectively. 



IV. EXPERIMENTAL RESULTS 

In the following we present our experimental results 
on the negative differential conductance phenomenon in 
gold-hydrogen nano junctions. The high stability atomic 
sized Au junctions were created by low-temperature me- 
chanically controllable break junction technique. The 
hydrogen molecules were directed to the junction from 
a high purity source through a capillary in the sample 
holder. The molecules were dosed by opening a solenoid 
needle valve with short voltage pulses, adding a typical 
amount of ~ 0.1 /imol of hydrogen molecules. 

The appearance of the NDC phenomenon cannot be 
generally attributed to definite parts of the conductance 
traces, but in AU-H2 junctions the NDC curves were 
quite frequently observed when the junction was closed 
after complete disconnection, as demonstrated by the 
trace in Fig. The conductance curve was measured 
by recording a large number of dI/dV(V) curves dur- 
ing a single opening - closing cycle, and extracting the 
differential conductance values both at zero bias and 
at Vdc = 120 mV. During the opening of the junc- 
tions no difference is observed, but during the closing 
of the junction the two curves show large deviation. The 
high bias trace resembles the traditional traces of pure 
gold junctions: during the approach of the electrodes 
an exponential-like growth is observed, but already at a 
small conductance value (< 0.1 Go) the junctions jumps 
to a direct contact with G = 1 Go- At low bias voltage 
the conductance grows to a much higher value (« 0.5 Go) 
before the jump to direct contact. The behavior of the 
conductance trace in Fig. 3^, agrees with the general trend 
shown by the conductance histograms (Fig. [3J)): at low 
bias a large variety of configurations is observed at any 
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FIG. 5: Panel (a) shows an experimental NDC curve (thick 
orange curve in the background) and the fit with a Dirac- 
delta distribution (open circles) and a finite width uniform 
distribution (thin black line). The fitting parameters are: 
E = 66.1 meV, N = 42, a = 0.15 Go, (Too = 0.0022 Go, : 
and T = 25/4.2 if, W = 0/4.5 meV for the Dirac-delta and 
the uniform distribution, respectively. Panel (b) shows exper- 
imental curves exhibiting smaller features between the NDC 
peaks. 



conductance value, whereas at high bias the weight in 
the region G = 0.1 — 0.8 Go is very small. Note that 
only a part of the configurations (a few percent of all 
traces) show NDC features, for other configurations irre- 
versible jumps are observed in the I— V curve at a certain 
threshold voltage. We have not found any evidence that 
the gold-hydrogen chains reported in our previous work 7 
would show NDC phenomenon or any peaklike structures 
in the differential conductance curve. 

Figure [1] has already presented an example of the huge 
negative differential conductance phenomenon. Using 
the finite temperature results of the asymmetrically cou- 
pled TLS model with a Dirac-delta distribution of the 
levels (W = 0) we have fitted our experimental data in 
Fig. [1] The fitting parameters are the energy, the tem- 
perature and the asymmetry parameter, whereas do and 
(Too are directly read from the experimental curve. As 
shown by the black solid lines in the figure, the model 
provides perfect fit to the data. The striking result of 
the fitting is the extremely large asymmetry parameter, 
N = 60. 

Figure shows an other example for the NDC phe- 
nomenon. For this curve the fitting with a Dirac-delta 
distribution provides a nonrealistic temperature of T — 
25 K. By fitting with a finite width uniform distribution 14 
of the excitation levels the temperature can be kept at the 
experimental value (4.2 K) resulting in a standard devi- 
ation of the upper states of W — 4.5 meV. Fig. shows 
further experimental curves exhibiting smaller features 
between the negative differential conductance peaks. In 
some measurements additional peaks appeared beside the 
main NDC peaks, which is demonstrated by the lower 
curve in the figure showing two smaller peaks at E = 37.6 
and 42.4 meV. These observations could be explained by 
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FIG. 6: The lower panel shows the variation on a NDC 
curve as a function of electrode separation. The curves were 
recorded during the closing of a disconnected junction. The 
upper panels show the variation of the fitting parameters as a 
function of the zero-bias conductance. The curves with zero 
bias conductances of 0.25, 0.35 and 0.45 Go correspond to the 
pushing of the electrodes by ~ 0.2A, ~ 0AA and « 0.6^4 
with respect to the initial curve with zero bias conductance of 
0.15 Go. The inset in the lower panel shows the distribution 
of the peak positions for 60 independent NDC curves. 



introducing further energy levels with different degener- 
acy, but such a model would be too much complicated 
without aiding the understanding of the general phe- 
nomenon. The upper curve shows a steplike decrease of 
~ 10% at E = ±39 meV resembling the vibrational signal 
of molecular junctions, although its amplitude is rather 
large for a vibrational spectrum. (The typical conduc- 
tance change due to vibrational excitations is 1 — 3% at 
the conductance quantum, 4 and it should vanish towards 
0.5 Go, where the negative point-contact spectroscopy 
signal turns to positive inelastic electron tunneling spec- 
troscopy signa h 12 ' 13 ) The observed steps could also be 
explained by the scattering on a symmetrically coupled 
TLS, for which the step-size can be as large as 50% of 
the zero bias conductance. 

The bottom panel in Fig. [6] shows the variation of the 
negative differential conductance curves by changing the 
electrode separation. The curves were recorded during 
the closing of a disconnected junction, similarly to the 
trace in Fig. [4^. The upper panels present the change of 
the fitting parameters as a function of the zero bias con- 
ductance. Both the excitation energy, E and the width 
of the excited levels, W decrease with increasing con- 
ductance, whereas the asymmetry parameter, N varies 
around a constant value. 

We have repeated our measurements for a large 
amount of curves showing NDC feature. The observed 
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FIG. 7: An illustration for the proposed model: a molecule 
strongly bound between the electrodes is excited to a loosely 
bound configuration, where different configurations have have 
similar binding energy. 



excitation energies have shown a broad distribution in 
the range 30 — 100 meV, as demonstrated by the inset in 
tho bottom panel of Fig. [5] Below 30meV we have not 
observed any peaks, whereas above 100 meV the junc- 
tions frequently become unstable, and the recording of 
reproducible I — V curves is not possible. In every case 
when the electrode separation dependence was studied E 
was decreasing by increasing zero-bias conductance show- 
ing a shift of the excitation energy by even a few tens of 
millivolts. The standard deviation of the excited levels 
was typically in the range 0.5 — lOmeV, also showing a 
decreasing tendency with increasing conductance. The 
asymmetry parameter was in the range rs 40 — 200. Note 
that with these parameters the average spacing of the ex- 
citation energies is smaller than the temperature smear- 
ing, thus the assumption for the continuous distribution 
of the levels is correct. 

In our experiments the NDC curves on gold-hydrogen 
junctions have not shown any telegraph noise within the 
bandwidth of our setup (100 kHz). 



V. DISCUSSION 

Our analysis of the two-level system models shows that 
in the case of two single levels the population of the up- 
per level remains very small in the narrow neighborhood 
of the excitation energy, which inhibits the appearance 
of sharp peaklike structures in the differential conduc- 
tance curve. By introducing a large asymmetry in the 
coupling constants, already a small voltage bias above 
the excitation energy causes a sudden flip of the occupa- 
tion numbers, resulting in peaks or even huge negative 
differential conductance in the dl/dV curves. 

The analysis of our experimental curves shows that the 
negative differential conductance phenomenon is success- 
fully fitted by the asymmetrically coupled TLS model, 
yielding extremely large asymmetry parameters (N m 
40 — 200). This result implies, that the molecular con- 
tact has a ground state with a well-defined molecular 
configuration, from which it can be excited to a large 
number of energetically similar excited states. A trivial 
explanation would be the desorption of a bound molecule 
to the vacuum, however it is hard to imagine, that after 
complete desorption a molecule always returns to exactly 
the same bound-state. As a more realistic explanation, 
a strongly bound molecule is excited into a large number 



of energetically similar loosely bound states, from which 
it can relax to the original bound state. 

A possible illustration is presented in Fig. \7\ In the 
ground state the molecule is strongly bound between the 
electrodes, whereas the excited states are loosely bound 
configurations where e.g. different angles of the molecule 
with respect to the contact axis have similar binding en- 
ergy, and also the molecule can diffuse to different sites 
at the side of the junction. This picture is consistent with 
the conductance trace in Fig. 2) at low bias the molec- 
ular contact is dominating the conductance, whereas at 
high bias the molecule is kicked out to the side of the 
contact, and the conductance trace resembles the behav- 
ior of pure gold junctions. It is noted though, that the 
pcaklike structures are observed with a large diversity of 
conditions, so in general Fig. [7] is regarded only as an 
illustrative picture 

The conclusive message of the analysis is, that the des- 
orption of a strongly bound molecule to a large number 
of energetically similar loosely bound states causes peak- 
like structures or even huge negative differential conduc- 
tance in the dl/dV curves. From the fitting of the ex- 
perimental curves only the large asymmetry parameter 
can be deduced, but the precise microscopic origin of 
the asymmetry cannot be determined without detailed 
microscopic calculations. Possible candidates are the dif- 
ferent arrangements and rotational states of the molecule 
with respect to the contact, but it is also very much prob- 
able, that the bound molecule is desorbed to different 
positions on the side of the contact, and then it can even 
diffuse away on the contact surface. In the latter case not 
necessarily the same molecule relaxes back to the initial 
bound state. All the above processes cause an effective 
asymmetry in the coupling. 

Similar results were obtained on tunnel junctions with 
higher resistance in Ref. 9|, where a hydrogen-covered Cu 
surface is studied with an STM tip. The authors present 
a phenomenological two-state model, which gives perfect 
fit to the data, but the microscopic parameters are hid- 
den in the model, and for the more detailed understand- 
ing of the results further analysis is required. The fitting 
with our asymmetrically coupled TLS model shows that 
the experimental curves in Ref. correspond to similarly 
large asymmetry parameters. In the case of an STM ge- 
ometry a direct picture can be associated with the asym- 
metry: a molecule bound between the surface and the tip 
can be excited to several equivalent states on other sites 
of the surface, away from the tip. The observation of the 
voltage dependent telegraph noise in Ref. [9| provides a 
direct measure of the occupation numbers, proving the 
population inversion in the system. 

The group of Jan van Ruitenbeek has demonstrated 
peaklike structures concentrating on contacts close to the 
conductance unit 4^ In this conductance range the rel- 
ative amplitudes of the peaks are much smaller than our 
NDC curves in the lower conductance regime, but the 
overall shape of the curves is similar. In some cases a 
direct transition from a steplike vibrational signal to a 
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FIG. 8: In the inset the vibrationally mediated TLS model is 
illustrated. The thick curve in the background is a fit to one 
of the experimental curves in Re], [g using the vibrationally 
mediated TLS model (the parameters are: o~o = 0.63 Go, <Ti = 
1.25 Go, huj = 40 meV, A = 3 meV, T = 7 K). The thin black 
line shows a fit with asymmetrically coupled TLS model with 
E = 41meV, T = 4.2 K, TV = 120 and W = 2.6 meV. 

peaklike structure was observed, which implies the cou- 
pling of the phenomenon to molecular vibrations. The 
coupling to the vibrational modes was also indicated by 
the isotope shift of the peak positions. A TLS model il- 
lustrated in the inset of Fig. [8] was proposed, in which the 
two base levels of the double well potential (with a split- 
ting of Aq) are separated by a wide potential well, which 
cannot be crossed. In both wells, however, a vibrational 
mode of the molecule can be excited (with E — Tilu), and 
at the vibrational level the potential barrier can already 
be crossed, causing a hybridization of the excited lev- 
els between the two wells. In this model the upper base 
level is unoccupied for eV < htv, but above the excitation 
energy the transition between the two wells is possible, 
and for ?wj A the upper base level suddenly becomes 
almost half occupied above the excitation energy. Simi- 
larly to the asymmetrically coupled TLS model, this vi- 
brationally mediated TLS model produces sharp peaklike 
structures in the dl/dV due to the sudden change of the 
occupation numbers. The authors claim that the cou- 
pling of a vibrational mode to a TLS magnifies the oth- 
erwise tiny vibrational signal, and thus the sharp peaks 
directly show the vibrational modes of the junction. 

Our analysis shows that the asymmetrically coupled 
TLS model proposed by us gives almost identical fits to 
the experimental curves in Ref. § as the vibrationally me- 
diated TLS model, which means that the automatic iden- 
tification of the peaklike structures with the vibrational 
modes of the junction is only possible if their growth from 
a vibrational spectrum is detected, and so the vibrational 
signal is anyhow resolved. In contrast, the negative dif- 
ferential curves demonstrated in our manuscript cannot 
be fitted with the model in Ref. |8| due to the large ratio of 
co/coo (~ 5 — 50). The authors use a symmetric coupling 
of the two levels, thus <7o/ cr oo < 2 for moderate voltages 
(eV > huj), and cto/coo < 3 for high bias (eV ^> TwS), 
where the third level starts to be populated as well. The 
model could be generalized by inserting asymmetric cou- 
pling, however in this case the inclusion of the vibrational 
level is not necessary any more for fitting the curves. The 
broad distribution of the NDC peak positions also indi- 



cates that the phenomenon is not related to vibrational 
energies. 

In our view the vibrationally mediated TLS model re- 
quires unique circumstances: the molecule needs to be 
really incorporated in the junction with well-defined vi- 
brational energies. Furthermore, two similar configura- 
tions are required, which have the same vibrational en- 
ergy, and for which the transition is only possible in the 
excited state. This phenomenon may be general in cer- 
tain well-defined molecular junctions, for which the tran- 
sition of a vibrational-like signal to a peaklike structure 
and the isotope shift of the peaks provides support. In 
contrast, the asymmetrically coupled TLS model can be 
imagined under a much broader range of experimental 
conditions, just a bound molecular configuration and a 
large number of energetically similar loosely bound (e.g. 
desorbed) states are required with some difference in the 
conductance. The fine details of the underlying physical 
phenomena may differ from system to system, however 
both models provide an illustration of a physical process 
that may lead to the appearance of peaklike structures 
in dl/dV curves of molecular junctions. 



VI. CONCLUSIONS 

In conclusion, we have shown that gold-hydrogen nano- 
j unctions exhibit huge negative differential conductance 
in the conductance range of 0.1 — 0.7 Go- The position of 
the peaks shows a broad distribution in the energy range 
30 — 100 meV. Similar features were observed in tunnel 
junctions with higher resistance using an STM setup, 9 
whereas in break-junctions with G ~ 1 Go peaklike struc- 
tures with smaller amplitude were detected^ These re- 
sults show that sharp peaklike structures are generally 
observed in the differential conductance curves of molec- 
ular nanojunctions under a wide range of experimental 
conditions: using different molecules, different electrode 
material and different contact sizes. All these experimen- 
tal results imply the explanation of the observations in 
terms of two-level system models. 

We have shown that a simple two-level system model 
cannot produce peaklike structures, but a TLS model 
with asymmetric coupling successfully fits all the experi- 
mental data. Our analysis shows that the peaklike curves 
correspond to large asymmetry parameters, which im- 
plies that a molecule from a bound state is excited to 
a large number of energetically similar loosely bound 
states. This picture provides a physical phenomenon that 
can appear under a wide range of experimental condi- 
tions and may explain the frequent occurrence of peaklike 
structures in dl/dV curves of various molecular systems. 

Our analysis also shows that a recently proposed vi- 
brationally mediated two-level system mode&i^ may be 
applicable for certain well-defined molecular configura- 
tions, but the general relation of the conductance peaks 
to vibrational modes is not possible without further ex- 
perimental evidences. 
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